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Abstract 

Nuclear matter and neutron stars are studied in the framework of an extended 
relativistic mean-field (RMF) model with higher-order derivative and density 
dependent couplings of nucleons to the meson fields. The derivative couplings 
lead to an energy dependence of the scalar and vector self-energies of the nu¬ 
cleons. It can be adjusted to be consistent with experimental results for the 
optical potential in nuclear matter. Several parametrisations, which give iden¬ 
tical predictions for the saturation properties of nuclear matter, are presented 
for different forms of the derivative coupling functions. The stellar structure 
of spherical, non-rotating stars is calculated for these new equations of state 
(EoS). A substantial softening of the EoS and a reduction of the maximum 
mass of neutron stars is found if the optical potential constraint is satisfied. 

Keywords: Relativistic mean-field model, Equation of state, Neutron stars. 
Density-dependent coupling, Derivative coupling. Optical potential 


1. Introduction 

The recent observation of two pulsars with approximately two solar masses 
[il,i presents a severe challenge to the theoretical description of cold high- 
density matter in /3-equilibrium. The equation of state (EoS) has to be suffi¬ 
ciently stiff in order to support such high masses of compact stars. Many models 
that are solely based on nucleonic (neutrons and protons) and leptonic (elec¬ 
trons and muons) degrees of freedom are able to reproduce maximum neutron 
star masses above two solar masses if the effective interaction between the nu¬ 
cleons becomes strongly repulsive at high baryon densities. However, additional 
hadronic particle species can appear at densities above two or three times the 
nuclear saturation density Usat ~ 0.16 fm“^. In most cases, these additional de¬ 
grees of freedom lead to a substantial softening of the EoS resulting in a reduced 
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maximum mass of the compact star below the observed values. This feature is 
well-known for models with hyperons - the so-called ’’hyperon puzzle”, see, e.g., 
ii and references therein - but was also observed in approaches that take 
excited states of the nucleons such as A(1232) resonances into account, see, e.g., 
@ and references therein. Usually, only specifically designed interactions can 
avoid the problem of too low maximum masses. 

Successful models of the baryonic contribution to the stellar EoS should be 
scrutinized whether they comply with other experimental constraints, e.g. with 
respect to the employed interactions. In the center of compact stars very high 
baryon densities are reached exceeding several times Usat and the corresponding 
Fermi momenta of the particles are much larger than those at saturation. This 
is particularly significant for models with only nucleonic degrees of freedom. 
Hence, not only the density dependence of the effective in-medium interaction 
between nucleons but also their momentum dependence becomes relevant. For 
densities near nsat this information is contained in the optical potential of nu¬ 
cleons that can be extracted from the systematics of elastic proton scattering 
on nuclei, see, e.g., 0, 0. A saturation of the real part of the optical potential 
is observed at high kinetic energies approaching 1 GeV. The momentum depen¬ 
dence of the in-medium interaction is also crucial in simulations of heavy-ion 
collisions Q. 

Typical approaches for the baryonic contribution to the stellar EoS are en¬ 
ergy density functionals that originate from nonrelativistic or relativistic mean- 
held models of nuclear matter. The most prominent cases among the former 
class are Skyrme energy density functionals, see reference 0 for an overview 
of different parametrizations. They are derived originally from the zero-range 
Skyrme interaction with a two-body contribution, which is an expansion up 
to second order in the particle momenta, and a density dependent three-body 
contribution, which is included in order to reproduce saturation properties of 
nuclear matter. Obviously, an extrapolation of the model to high momenta is 
questionable given the limited form of the momentum dependence. Examples 
of the latter class, frequently denominated covariant density functionals, can be 
inferred from relativistic Lagrangian densities. In conventional models, see ref¬ 
erence ll| for a wide collection of different parametrizations, a nucleon optical 


potential in the medium can be derived from the relativistic scalar and vector 
self-energies, see section [S] It exhibits a linear increase with energy, which is 
in contradiction with the expectation from experiment. In general, one would 
expect that the nucleon self-energies itself depend explicitly on the particle mo¬ 
mentum or energy as, e.g., in Dirac-Brueckner calculations of nuclear matter 


12|. However, this is not realized in standard RMF models. 

There are particular extensions of relativistic mean-field (RMF) models that 
contain nucleon self-energies with an explicit energy or momentum dependence. 
This dependence cannot be introduced in a relativistic model in a simple para¬ 
metric form because it affects, e.g., the definition of the conserved currents. In 
extended systematic approaches new derivative couplings between the nucleon 
and meson fields are introduced that allow to reproduce the energy dependence 
of the optical potential as extracted from experiments. One of the earliest 
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RMF models with scalar derivative couplings was presented in reference [13| . A 
rescaling of the nucleon fields removed the explicit momentum dependence of 
the self-energies but lead to a considerable softening of the EoS. More general 
couplings of the mesons to linear derivatives of the nucleon fields were considered 
in reference 1^ with an application to uniform nuclear matter. With appropri¬ 
ately chosen coupling constants a reduction of the optical potential was found as 
compared to the strong linear energy dependence in conventional RMF models. 
The model was further extended in reference [l5| assuming a density depen¬ 
dence of the couplings. It was successfully applied to the description of finite 
nuclei. Notable new features were the increase of the effective nucleon masses 
(usually rather small in order to explain the strong spin-orbit interaction in 
nuclei) and correspondingly higher level densities close to the Fermi energy in 
nuclei in better concordance with expectations from experiments. Though, us¬ 
ing couplings only linear in the derivatives leads to a quadratic dependence of 
the optical potential on the kinetic energy with a decrease for energies exceeding 
1 GeV. Couplings to all orders in the derivative of the nucleons were introduced 
in the so-called nonlinear derivative (NLD) model [lBj[l3| assuming a particular 
exponential dependence on the derivatives but no density dependence of the 
nucleon-meson couplings or self-couplings of the mesons. The general formal¬ 
ism was developed and applied to infinite isospin symmetric and asymmetric 
nuclear matter with a particular choice of the non-linear derivative terms that 
lead to an energy dependence of the self-energies. In reference [I^ the approach 
was extended with generalized non-linear derivative couplings of any functional 
form in the field-theoretical formalism allowing for a momentum or energy de¬ 
pendence. Nonlinear self-couplings of the a meson field were added in order 
to improve the description of characteristic nuclear matter parameters at sat¬ 
uration. The application of this version of the NLD model to stellar matter 
yielded a maximum neutron star mass of 2.03 Mgoi barely satisfying the obser¬ 
vational constraints but the dependence of the result on the model parameters 
was not explored in detail. The NLD model was also applied to the description 
of bulk properties of nuclear matter in reference [l^. A softening of the EoS 
and maximum masses of neutron stars substantially below 2 Mgoi were found 
with different parametrizations assuming an energy dependence of the couplings 
but nonlinear self-couplings of the mesons or density dependent meson-nucleon 
couplings were not considered. Properties of finite nuclei were studied in refer¬ 
ence [20| after adding meson self-interactions in the Lagrangian. A qualitative 
description similar to conventional RMF models was achieved but neutron star 
properties were not examined in this extended model. 

In this work we introduce a more flexible extension of the nonlinear deriva¬ 
tive model assuming density dependent meson-nucleon couplings in addition. 
Instead of using derivative operators that generate an explicit momentum de¬ 
pendence of the self-energies, we will use a functional form that leads to an 
energy dependence. This approach will also be more suitable for a future ap¬ 
plications of the DD-NLD approach to nuclei since the relevant equations and 
their numerical implementation are simplified. Here, the equations of state of 
symmetric and asymmetric nuclear matter will be calculated for different choices 
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of the derivative coupling operators that lead to a saturation of the optical po¬ 
tential at high energies as derived from experiments. They are compared to 
the results of a standard RMF model with density dependent couplings that is 
consistent with essentially all modern constraints for the characteristic nuclear 
matter parameters at saturation. The parameters of the DD-NLD models are 
chosen such that these saturation properties are reproduced. The effect of the 
optical potential constraint on the mass-radius relations of neutron stars will be 
studied. 

The paper is organized as follows: In section[5]the Lagrangian density of the 
DD-NLD approach is presented. The field equations in mean-field approxima¬ 
tion and the energy-momentum tensor will be derived. The relevant equations 
for the case of infinite nuclear matter will be considered in more detail in section 
[3l The parametrization of the density dependent couplings and the functional 
form of the derivative coupling functions is discussed in sectional Results for the 
energy dependence of the optical potential, the EoS of nuclear matter and the 
mass-radius relation are presented in section [S] for various versions of the model. 
Conclusions are given in section [51 Detailed expression for various densities are 
collected in [Appendix A| 

2. Lagrangian density and field equations of the DD-NLD model 

In most RMF models the effective interaction between nucleons is described 
by an exchange of mesons. Usually, cr and w mesons are introduced to consider 
the attractive and repulsive contributions to the nucleon-nucleon potential, re¬ 
spectively. They are represented by isoscalar Lorentz scalar and Lorentz vector 
fields a and In order to model the isospin dependence of the interaction, the 
exchange of p mesons is included. It is denoted by the isovector Lorentz vector 


field pfi in the following. The Lagrangian density in the DD-NLD approach^ 


( 1 ) 


£ = £, 


+ £] 


+ A 


'int 


nuc 


•mes 


contains contributions of the free nucleons = (4'p, 4'„) with mass m 



( 2 ) 


in a symmetrized form and of free mesons 



( 3 ) 


with the field tensors 


= 9tJ.uJu - and - di,p^ 


( 4 ) 


^Natural units with h = c = 1 are used in the following. 
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of the isoscalar w meson and the isovector p meson, respectively. The arrows in 
equation ([2]) denote the direction of differentiation. 

Standard RMF models assume a minimal coupling of the nucleons to the 
meson fields leading to 

£i„t = (5) 

for the interaction contribution to the total Lagrangian density C with meson- 
nucleon couplings Ti (i = a, Lu, p). We assume that they depend on the vector 
density riy, see equation (I25p for the explicit definition. In the derivative cou¬ 
pling model the nucleon field (^') is replaced in £int by with op¬ 

erator functions which can be different for the various mesons m = a,ui,p. 
They can be expanded in a series 

OO 7(^) 

Vm{x) = ^-^X^ ( 6 ) 

n\ 

n—O 

with numerical coefficients . The argument x contains derivatives idp that 
act on the nucleon field. More specifically we write 

X = v^idp — sm (7) 

as a hermitian Lorentz scalar operator with an auxiliary Lorentz vector = 
(vo,v) and a scalar factor s. Hence the interaction contribution in the NLD 
model is written as 

£i„t = + (8) 

in a symmetrized form with respect to the derivative operators i-e., 




(9) 


k^O 


i — 
T’m 

OO 

/c=0 

(10) 


with coefficients 



Obviously, no derivatives appear for the choice T>m = 1 (corresponding to = 
d„o) and the standard form ([5]) is recovered. 
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When spatially inhomogeneous systems with Coulomb interaction are con¬ 
sidered, £mes and £int can be complemented with the appropriate contributions. 
Since only uniform matter is considered in the following, we do not give them 
here explicitly. 

The field equations of nucleons are derived from the generalized Euler- 
Lagrange equation 


dL 

difr 


i=l 


i-yd^ 


dC 


d{da 


= 0 


iV’r 


( 12 ) 


for For the meson fields (pr = cr, the standard Euler-Lagrange 

equation applies, i.e. only the i = 1 term in equation m is relevant since 
higher-order derivatives of the meson fields do not appear in the Lagrangian 
density C. 

Details can found in references The Dirac equation 

[ 7 ^ {idf^ - - (m - E)] ^ = 0 (13) 


for the nucleons looks formally the same as in standard RMF approaches but 
the scalar (E) and vector (E^) self-energy operators now contain the derivative 
operators Dm- They are given by 

E = r„aV„ (14) 


and 








with the ’rearrangement’ contribution 

,1 


K = - 


n- 

+r; 

-r 


>7 ( 

(■ 

>5 (*^' 




(15) 

(16) 


'i’V p'jjyT'i’ + \E' 7 i 




4' 


containing derivatives 


r' = 


4' -f 4'D„4' 


dT,, 


dn„ 


(17) 


of the coupling functions. In the case of inhomogeneous systems and a non¬ 
vanishing three-vector component v of the auxiliary vector , additional contri¬ 
butions in 0 and (IT5|) will appear. In the present application of the DD-NLD 
model, however, we do not consider this case. The field equations of the mesons 
are found as 

1„ 


dpd^a + mla = -T^ 

(^^7^4'+¥7"^^4- 
ipp + 'I'T7 




2 u 

■ m,„u) = 


^^p4' 
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(18) 

(19) 

( 20 ) 





with source terms containing derivative operators. 

The conserved baryon current in the DD-NLD model is given by 

^ ( 21 ) 

i—p,n 


with the norm operator 

+ T^a (d^V,) - r^a;„7“ (d^^u.) - rpP„7“^ (22) 

where d^T>m is the derivative of operator with respect to the momentum 
Pp. = i.e. 

CO 

= ^ ( 23 ) 

k=l 

and (...) denotes the summation over all occupied states. The current ([^ is 
not identical to the vector current 

= E (^*7'^'!'*), (24) 

i—p^n 


which is used to define the vector density 



appearing as the argument of the coupling functions Ti. The energy-momentum 
tensor assumes the form 

rj^flU _ E (4-,iV>"4-,)-g^"(£). (26) 

i—p,n 


Then the energy density s and pressure p are found from e = and p = 
X)i=i 7"“/3j respectively. 

3. DD-NLD model for nuclear matter 

In the case of stationary nuclear matter, the equations simplify considerably 
since the system is homogeneous and the meson fields, which are treated as 
classical fields, are constant in space and time. Positive-energy solutions of 
the Dirac equation (USD are plane waves 'ii = Ui exp (—for protons and 
neutrons with Dirac spinors Ui, which are normalized according to 

= u,N°u, = 1 (27) 

with the time component of the norm operator (I22|) . They depend on the 
effective mass 

m* = rrii — Yji (28) 
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and effective momentum 


(29) 


related by the dispersion relation 

P^pI = «)" ■ (30) 

The derivative id^ in the Dm operators can be replaced by the corresponding 
four-momentum = {Ei,pi) resulting in a simple function Dm depending on 
the energy Ei and the momentum pi of the nucleon. 

Using the identity 


)] (31) 


the conserved current and the energy-momentum tensor can be written as 


i—p,n 


<Pp nf 

(27r)3 n« 


(32) 


and 


i—p,n 


(Pp 




respectively, with the four-momentum 





(33) 


(34) 


and spin degeneracy factors Ki = 2. The integration runs over all momenta p 
with modulus lower than the Fermi momenta ppi in the no-sea approximation. 
They are defined through the individual nucleon densities 




(35) 


Without the preference for a particular direction in infinite nuclear matter, 
the spatial components of the Lorentz vector meson fields vanish and the aux¬ 
iliary vector in equation ([7|) is set to = 5po such that the Dm functions only 
depend on the nucleon energy Ei. Without isospin changing processes, only the 
third component of the isovector p field has to be considered in the field equa¬ 
tions for the mesons. Using the abbreviations uj = uP and p = the meson 
fields are immediately obtained from 


a 


LO 


P 


^ i—p,n 




m- 

r„ 


mf: 


i—p.n 


r 


— Up = 

mi m 




P 2 = 


i—p,n 


(36) 

(37) 

(38) 







with source densities ria, ria,, and Up. The self-energies simplify to 


Si = T^aDo- (39) 

Si = ^pT^,ipDp -\- S^ (40) 

s, = 0 (41) 

with T'ip = 1 (—1) for protons (neutrons) and the ’rearrangement’ contribution 

S^j = + T'ppup - r'^aua (42) 

is independent of the nucleon energy. The dispersion relation reads 

El = \J+ {rrii - Sif + Vi (43) 

if we introduce the energy-dependent scalar potentials Si{E) = Si and vec¬ 


tor potentials Vi{E) = S-*. Explicit expressions for the various densities and 
thermodynamic quantities of the DD-NLD model are given in [Appendix A 

4. Parametrization of the DD-NLD model 

For the application of the DD-NLD model to nuclear matter the parameters 
need to be specified. Besides the usual parameters of a RMF model with density 
dependent couplings the form of the Dm functions has to be given. We assume 
identical functions for all mesons, i.e. D = D,j = = Dp, and consider three 

functional dependencies: 

D1 a constant D = 1, which corresponds to a usual RMF model with density 
dependent couplings, 

D2 a Lorentzian form D = 1/(1 + x'^), 

D3 an exponential dependence D = exp {—x) 

with X = {Ei — mi)/A because we set = 6po and s = 1 in equation ([7]). The 
parameter A regulates the strength of the energy dependence. For the proton 
and neutron masses the experimental values of nip = 938.272046 MeV/c^ and 
nin = 939.565379 MeV/c^, respectively, are used. The meson masses are set to 
TOcr = 550 MeV/c^, rriuj = 783 MeV/c^, and nip = 763 MeV/c^. The density 
dependence of the meson nucleon couplings has the same form as introduced 
in reference [2lj. We assume a dependence on the vector density Uy as defined 
in equation (1251) , which is not identical to the zero-component of the conserved 
baryon current (EH- See [Appendix A| for explicit expressions. This choice 
simplifies the rearrangement terms considerably. 

For the isoscalar mesons m = a, ui the coupling is written as 

fm{x) (^4) 
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Table 1: Parameters of the meson coupling functions for three choices of the D functions and 
different values of cut-off parameter A. 


meson 

parameter 

D1 

D2 


D3 


A [MeV] 

- 

400 

500 

600 

700 

a 

ro-('^ref) 

10.72913 

10.93466 

10.86315 

9.74679 

9.89158 


(Xfj 

1.36402 

1.35816 

1.36015 

1.38410 

1.38064 


ba 

0.53404 

0.51914 

0.52433 

0.61515 

0.60127 


C-a 

0.86714 

0.83989 

0.84931 

1.00615 

0.98211 


dfj 

0.62000 

0.62998 

0.62648 

0.57558 

0.58258 

LU 

rLj('^ref) 

13.29858 

13.56462 

13.47215 

12.0503 

12.23457 


(X(jj 

1.3822 

1.3822 

1.3822 

1.3822 

1.3822 



0.42253 

0.42253 

0.42253 

0.42253 

0.42253 


Cuj 

0.71932 

0.71932 

0.71932 

0.71932 

0.71932 


duj 

0.60473 

0.68073 

0.68073 

0.68073 

0.68073 

p 

rp(^ref) 

3.59367 

3.67852 

3.64957 

3.18819 

3.25233 


Qp 

0.48762 

0.48954 

0.48872 

0.34777 

0.36279 


Tlref [fm“3] 

0.15000 

0.14618 

0.147485 

0.16515 

0.16268 


with functions 


fm{^) 


1 H“ 6777,(x H- dm) 

1 Cm{^ dm)‘^ 


(45) 


that depend on the argument x = nyjriret and contain coefficients Cm, 

and dm- For the p-meson coupling we set 


rp(n„) = rp(nref) exp [-Op {x - 1)] . (46) 

In order to reduce the number of independent parameters we demand that the 
conditions /cr(l) = /tj(l) = 1 and /"(O) = /"(O) = 0 hold. Hence, there are only 
two independent coefficients in the functions fm for each of the isoscalar mesons. 
The overall magnitude of the couplings is given by the couplings rm(nref) at a 
reference density Uref- We require that the characteristic saturation properties 
for the three choices of the D function are identical and close to current values 
extracted from experiments. In particular, we set the saturation density to 
nsat = 0.15 fm“^, the binding energy per nucleon at saturation to H = 16 MeV, 
the compressibility to K = 240 MeV, the symmetry energy to J = 32 MeV 
and the symmetry energy slope coefficient to L = 60 MeV. Furthermore, we set 
the effective nucleon mass at saturation to mes = 0.5625 Wnuc (related to the 
strength of the spin-orbit potential in nuclei) and fix the ratios = 

—0.15 and /"(1)//(),(1) = —1.0 in order to determine the coefficients in the 
functions fm uniquely. These values are close to those of the parametrization 
DD2 [ 2 ^ that was fitted to properties of nuclei and predicts a neutron star 
maximum mass of 2.4 Mgoi- 
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Figure 1: The optical potential C/opt as a function of the kinetic energy -Fkin = E — rUnuc of a 
nucleon in symmetric nuclear matter at saturation density in RMF models with parametriza- 
tions D1 and D2 compared to two fits from Dirac phenomenology. See text for details. 


Explicit values of the model parameters are given in tabled] with two choices 
of the cut-off parameter A for the cases of Lorentzian and exponential func¬ 
tions D. Note that the coefficients of the function are identical for all five 
parametrizations due to the constraints. The reference density Uref is not nec¬ 
essarily identical to the saturation density risat because in the case of explicit 
derivative couplings the vector density Uy is different from the conserved baryon 
density ub = = np + Un- 


5. Results 


The nonlinear derivative couplings are introduced in the RMF model in order 
to improve the energy dependence of the optical potential Uopt- The elastic 
proton scattering on nuclei of different mass number A can be well described in 
Dirac phenomenology with scalar (S) and vector {V) potentials, which smoothly 
vary with A and the energy of the projectile [71, l8|. From these global fits the 
optical potential in symmetric nuclear matter at saturation density is obtained 
as a function of the kinetic energy Fkin = E — mnuc in the limit A ^ oo. There 
are different definitions of the nonrelativistic optical potential when it is derived 
from relativistic scalar and vector self-energies. Here we use the form 


p c2 _ ^2 

Uopt{E) = - 1 / - 5 + ^ 

TRnuc -^^nuc 


(47) 


with S = Ep and F = as in references 14, 1^ 16, 17, 1^. In conventional 
RMF models without derivative couplings, the scalar and vector potentials are 
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Figure 2: The optical potential C/opt as a function of the kinetic energy -Ekin of a, nucleon 
in symmetric nuclear matter at saturation density in RMF models with parametrizations D1 
and D3 compared to two fits from Dirac phenomenology. See text for details. 


constant in energy and the optical potential (j47p is just a linear function in 
energy. This is clearly seen in figures [I] and [5] as a full black line for the calcu¬ 
lation with the parametrization Dl. In contrast, the optical potentials derived 
from the scalar and vector potentials in Dirac phenomenology from two differ¬ 
ent fits 0 are much smaller at high energies and exhibit a saturation for i?kin 
approaching 1 GeV. At low energies, the optical potentials from experiment 
behave more similar as that of the theoretical model concerning the absolute 
strength and the energy dependence. In figure [T] ([2]) the result for the DD-NLD 
parametrization D2 (D3) is depicted for two values of the cut-off parameter A. 
Here, a reasonable description of the experimental optical potential is achieved 
due to the energy dependence of the nucleon self-energies. The difference be¬ 
tween parametrizations D2 and D3 is not very significant. The dependence of 
Uopt on A is stronger for D2. The deflection of the DD-NLD curve for that 
of the standard RMF model Dl appears at lower kinetic energies for the D3 
parametrization as compared to the D2 case. In the DD-NLD model, the op¬ 
tical potential can be calculated easily for other baryon densities and arbitrary 
neutron-proton asymmetries. Because there are no experimental data available 
for these general cases, we refrain from presenting the results here. But see 


14j for the systematics with density in the linear derivative coupling 


reference 
model. 

The reduction of the optical potential at high kinetic energies, which orig¬ 
inates from the energy dependence of the self-energies, is also reflected in the 
equation of state. In figure|3]the energy per nucleon E/A (without the rest mass 
contribution) is depicted as a function of the baryon density ub in symmetric 
nuclear matter (left panel) and neutron matter (right panel). In both cases. 
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rig [fm'^] iig [fm'^] 


Figure 3: Energy per nucleon EjA as a function of the baryon density ub for symmetric 
nuclear matter (a) and pure neutron matter (b). Results are given for three different choices 
of the D function and different values of the cut-off parameter A. 


a substantial softening of the EoS is found as compared to the standard RMF 
calculation with parametrisation Dl. The effect is stronger for an exponential 
energy dependence of the self-energies (D3) than for the case of a Lorentzian de¬ 
pendence (D2). By construction, all EoS are identical at the saturation density 

ttsat ■ 

The DD-NLD model can be used to predict the properties of neutron stars. 
Here, the EoS of stellar matter is required. It is obtained by adding the con¬ 
tribution of electrons to the energy density and pressure of the baryons. The 
conditions of charge neutrality and /3 equilibrium fix the lepton density and 
proton-neutron asymmetry uniquely. Since the present model calculations treat 
only homogeneous matter, a suitable EoS for the crust of neutron stars has to be 
added at low densities. We use the standard Baym-Pethick-Sutherland (BPS) 
crust EoS [ 2 ^. The mass-radius relation of neutron stars is found finally by 
solving the Tolman-Oppenheimer-Volkoff equations 2J,l25j. It is shown for the 
five models of this work in figure0]together with the masses of the two most mas¬ 
sive pulsars observed so far. The model without an energy dependence (black 
full line) can explain without any difficulties the large neutron star masses from 
astrophysical observations [l|,[2] because the EoS is rather stiff at high densities. 
In contrast, for the EoS of neutron star matter in the DD-NLD models that 
are consistent with the optical potential constraint, a serious reduction of the 
maximum neutron star mass is seen. These models have even problems to reach 
typical masses of about 1.4 Mgoi of ordinary neutron stars. Deviations from the 
predictions of the standard model Dl start to appear already at masses below 
0.7 Msoi. There is also an effect on the neutron star radius, which is found to 
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Figure 4: Mass-radius relation of neutron stars for different choices of the D functions and 
cut-off parameters A in the DD-NLD model. The two shaded bands refer to astrophysical 
mass measurements of the pulsars PSR J1614 — 2230 [3] and PSR J0348 -I- 0432 Q. 


be smaller in the D2 and in the D3 model. Explicit values for the maximum 
mass as well as the radius and central density at this extreme conditions are 
given in tabled For the DD-NLD models D2 and, in particular, D3 the central 
densities is a star of maximum mass are considerably higher than those of the 
standard RMF model without an energy dependence of the couplings. 

It is worthwhile to compare our results for the mass-radius relation of neutron 
stars with those of previous versions of the NLD model. In the approach of 
reference 18| it was possible to reach a maximum neutron star mass of about 
2 Msoi with their choice of parameters. In contrast, the results of reference 
[l9l| using different functional forms of the couplings indicate a reduction of the 
maximum mass below the observed values in line with our calculations. All 
three versions of the NLD model are adjusted to similar values of the nuclear 
matter parameters at saturation, such as saturation density, binding energy, 
compressibility or symmetry energy consistent with experimental constraints. 
Nevertheless, the predictions for matter properties at supra-saturation densities 
are rather different due to the various choices in the models to represent the 
effective in-medium interaction. In our paper, an explicit density dependence 
of the couplings is considered whereas in [l8j nonlinear self-couplings of the 
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Table 2: Maximum mass, corresponding radius and central density of neutrons stars in the 
DD-NLD models with different parametrizations. 


model 


Almax [AIsol] 

^(Mmax) [km] 

^central(-^inax) ] 

D1 


2.37 

11.57 

0.88 

D2, A = 

500 MeV 

1.48 

11.75 

0.95 

D2, A = 

400 MeV 

1.36 

11.63 

0.97 

D3, A = 

700 MeV 

1.26 

10.13 

1.41 

D3, A = 

600 MeV 

1.19 

9.99 

1.46 


meson fields were assumened. In the approach of [l^ to nuclear matter neither 
meson self-couplings nor a density dependence of the couplings were used. A 
reasonable description of nuclear matter near saturation does not determine the 
high-density behaviour of the EoS uniquely, in particular due to the additional 
freedom with the explicit momentum/energy dependence in the NLD approach 
as compared to standard RMF models. Of course, the maximum mass constraint 
could be used directly in the determination of the model parameters. But also 
additional constraints at high densities, e.g. from heavy-ion collisions, might 
help to reduce the uncertainties in the extrapolation from low to high densities 
in the future. 


6. Conclusions 

There are several aspects that have to be taken into account in constraining 
models of dense matter for the application to neutron stars. In phenomenolog¬ 
ical models, the characteristic saturation properties of nuclear matter and the 
density dependence of the effective interaction are usually addressed. However, 
less attention is paid to its energy or momentum dependence. 

Introducing non-linear derivative couplings into RMF models, it is possible 
to generate an energy dependence of the nucleon self-energies such that the 
optical potential in nuclear matter, which is extracted in Dirac phenomenology 
from elastic proton-nucleus scattering experiments, can be well described up to 
energies of 1 GeV. Considering density-dependent nucleon-meson couplings at 
the same time, a very flexible model is obtained. Its parameters can be fitted to 
the usual nuclear matter constraints even for different functional forms of the 
energy dependent couplings. 

In the current version of the DD-NLD model, the energy dependence of the 
self-energies causes a softening of the EoS at high densities, both for symmet¬ 
ric nuclear matter and pure neutron matter. This effect is independent of the 
appearance of additional degrees of freedom, such as hyperons or deltas. As a 
result, it becomes more difficult to obtain very massive neutron stars consistent 
with the observational constraints. In the present work, the density dependence 
of the w meson coupling was kept fixed in the parametrizations and only a few 
choices for the energy dependence were tested. Only constraints near the nu- 
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clear saturation density were used to determine the model parameters. There 
remains substantial freedom in the model to be explored in order to find a suit¬ 
able parametrization that is consistent with the optical potential and maximum 
neutron mass constraint. Nevertheless, the results of the our study indicate 
that the optical potential constraint has to be taken seriously into account in 
the development of realistic phenomenological models for dense matter. 

RMF models with self-energies that explicitly depend on the nucleon mo¬ 
mentum or energy can be applied to the simulation of heavy-ion collisions using 
relativistic transport approaches. Here the equation of state can be tested at 
supra-saturation densities. It is well known that an energy/momentum depen¬ 
dence of the effective in-medium interaction is mandatory for a proper descrip¬ 
tion and analysis of experimental data. Such an approach can help to constrain 
the parameters of the present model at densities that are not accessible in the 
description of finite nuclei. 

In our work, an explicit energy dependence of the nucleon-meson couplings 
was favored. It allows to apply the DD-NLD approach to the description of 
nuclei without major difficulties. Such an additional investigation of the model 
will permit a better control on the parameters. In particular, the interplay be¬ 
tween the choice of the functional form of the energy dependence, the cutoff 
parameters and the density dependence of the couplings can be studied with its 
impact for the prediction of maximum neutron star masses. With a larger num¬ 
ber of observables as constraints it can be expected that firmer conclusions can 
be drawn on the compatibility of nuclear matter and finite nuclei descriptions 
in the DD-NLD model. Work in this direction is in progress. 
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Appendix A. Densities and thermodynamic quantities in the DD- 


NLD model 


The proton and neutron densities ([M]) in the DD-NLD model are easily 
calculated through a momentum integration of the relevant integral. For other 
quantities, however, it is more convenient to introduce an energy integration 
by substitution because the scalar (Si) and vector (I^) potentials are explicit 
functions of the energy Ei of the nucleon i. With the dispersion relation (143 1) 
we obtain 



(A.1) 


and the derivative 



(A.2) 


16 





with 


dEi 


(r^w + Fpra^jp) 


dD 

dE,, 


and 


dEi dEi ■ 


Introducing the scalar source densities 


n. 


{sD) 


^(max) 


= 7r\ dE, p{Ei) [m, - 5,(£;,)] Z?(£;,) 


(A.3) 

(A.4) 


(A.5) 


and vector source densities 


(vD) 

n) 


rPFi j 3„ p * 

'‘•X 

„{max) 

^ / dA, p(A,) [A - , 


the total source densities 


TLa 

= 4^^)+ 4*^) 

(A.7) 

riu, 

= 

(A.8) 

Up 

- _ (vD) 

— ILp IL^ 

(A.9) 


in the field equations dMl), (Ezl), and (l38l) are found. The lower and upper 
boundaries of the integrals are determined by solving the equations 


^(min) _ 


+ix(4““0 


(min)N 


and 


E, 


(max) 


= \ PF^ + 




(max) 


(A.IO) 


(A.11) 


respectively, with the Fermi momenta ppi from equation ()35p . The argument 

(v) (v) 

'P “t“ Ce 

d^p E* 


(v) (v) 

of the coupling functions riy = rip + rin can be obtained from 


{v) 

n) = Ki 


(2^)3 no 


(A.12) 


^(max) 


27r3 / rCmin) 


dE, p{Ei) [A, - F,(A,)] . 
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The energy density assumes the form 


% -p,7fc 

^(max) 



(A.13) 


^ r 

2tT^ _/^(min) 


L 


dE,p{E,)U°iE,)E,-{C) 


and the pressure is given by 


P 



(A.14) 


with 



Using a partial integration and equation (jA.2|) . the thermodynamic identity 





(A.16) 


with the chemical potential pi = Ei{ppi) is easily confirmed. 
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